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Abstract 

As a benchmark exercise for testing software and methods developed in Hycon2 for decen- 
tralized and distributed control, we address the problem of designing the Automatic Genera- 
tion Control (AGC) layer in power network systems. In particular, we present three different 
scenarios and discuss performance levels that can be reached using Centralized Model Predic- 
tive Control (MPC). These results can be used as a milestone for comparing the performance 
of alternative control schemes. Matlab software for simulating the scenarios is also provided 
in an accompanying file. 



*The research leading to those results has received funding from the European Union Seventh Framework Pro- 
gramme [FP7/2007-2013] under grant agreement n° 257462 HYCON2 Network of excellence. 
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1 Introduction 



An example of a real application that can benefit of decentralized and distributed control schemes 
is the regulation of a Power Network System (PNS). We consider a PNS as composed by several 
power generation areas coupled through tie-lines [Saa02]. The aim is to design the Automatic 
Generation Control (AGC) layer for frequency control with the goal of: 

• keeping the frequency approximately at the nominal value; 

• controlling the tie-line powers in order to reduce power exchanges between areas. In the 
asymptotic regime each area should compensate for local load steps and produce the required 
power. 

We consider thermal power stations with single-stage turbines. The dynamics of an area 
equipped with primary control and linearized around equilibrium value for all variables can be 
described by the following continuous-time LTI model [Saa02] 



= A u x[ t ] + BiU[i] + L t AP Li 



■ 2^ A v x U] 



(1) 



where = (A0j, Aw*, AP mi , AP Vi ) is the state, um = AP re f i is the control input of each area, 
APl is the local power load and A/i is the sets of neighboring areas, i.e. areas directly connected 
to ££ through tie-lines. The matrices of system (1) are defined as 
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(2) 



For the meaning of constants as well as some typical parameter values we defer the reader to Table 
1. 

We note that model (1) is input decoupled since both AP re f i and APi Ji act only on subsystem 
E£ . Moreover, subsystems are parameter dependent since the local dynamics depends on the 

quantities ~Yh. ■ 

In the following we introduce three scenarios corresponding to different interconnection topolo- 
gies of generation areas. The model parameters and constraints on A0i and on AP re f i for systems 
in all Scenarios are given in Table 2. We highlight that all parameter values are within the range 
of those used in Chapter 12 of [Saa02]. We define M as the number of areas in the power network. 
For each scenario, discrete-time models £[,•] with T s = 1 sec sampling time are obtained from £j?j 
using two alternative discretization schemes. 



Exact discretization of the overall system (acronym D); 



APr 



• Discretization system-by-system, i.e. exact discretization for each area treating ui 
and xu], j € A/i as exogenous inputs (acronym Dss). 

In particular, we note that Dss preserves the input-decoupled structure of ES while D does not. 
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Deviation of the angular displacement of the rotor with respect to the stationary reference axis on the stator 
Speed deviation of rotating mass from nominal value 
Deviation of the mechanical power from nominal value (p.u.) 
Deviation of the steam valve position from nominal value (p.u.) 
Deviation of the reference set power from nominal value (p.u.) 
Deviation of the nonfrequency-sensitive load change from nominal value (p.u.) 

t , ■ , , „ , „ kinetic energy at rated speed /,-,,,. n mi \ 

Inertia constant denned as Hi = A 1 t-. - (typically values m range 1 — 10 sec) 

machine rating w 1 J ° 1 1 ' 

Speed regulation 

x-. n rl percent change in load 

ucniiccl as — \ = — ? 

change in frequency 

Prime mover time constant (typically values in range [0.2 — 2] sec ) 
Governor time constant (typically values in range [0.1 — 0.6] sec ) 
Slope of the power angle curve at the initial operating angle between area i and area j 



Tabic 1: Variables of a generation area with typical value ranges [Saa02]. (p.u.) stands for "per 
unit" . 





Area 1 


Area 2 


Area 3 


Area 4 


Area 5 


Hi 


12 


10 


8 


8 


10 


Ri 


0.05 


0.0625 


0.08 


0.08 


0.05 


A 


0.7 


0.9 


0.9 


0.7 


0.86 


T u 


0.65 


0.4 


0.3 


0.6 


0.8 


T gi 


0.1 


0.1 


0.1 


0.1 


0.15 





Area 1 


Area 2 


Area 3 


Area 4 


Area 5 




lk[l,i]||oo < 0.1 


||Z[2,1]II°° ^ 0- 1 


l|Z[3,l]H°° ^ O' 1 


||£[4,1]IU ^ 0- 1 


l|^[5,l] lloo < 0.1 




|m[i]||oo < 0.5 


I |W[2] I |oo < 0.65 


I |w[3] 1 loo < 0.65 


\\u[4] oo < 0.55 


||«[5]||oo < 0.5 



P\2 — 4 P23 — 2 P34 — 2 P45 — 3 P25 — 3 
Table 2: Model parameters and constraints for systems Em, i e 1, . . . , 5. 



1.1 Scenario 1 

We consider four areas interconnected as in Figure 1. We will simulate Scenario 1 using the load 
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Figure 1: Power network system of Scenario 1 
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steps specified in Table 3. 



Step time 


Area i 




5 


1 


+0.15 


15 


2 


-0.15 


20 


3 


+0.12 


40 


3 


-0.12 


40 


4 


+0.28 



Table 3: Load of power AP^ (p-u.) for simulation in Scenario 1. +AP^ means a step of required 
power, hence a decrease of the frequency deviation Awi and therefore an increase of the power 
reference AP r e/ ; . 



1.2 Scenario 2 



We consider the power network proposed in Scenario 1 and add a fifth area connected as in Figure 
2. We will simulate Scenario 2 using the load steps specified in Table 4. 
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Figure 2: Power network system of Scenario 2 



Step time 


Area i 




5 


1 


+0.10 


15 


2 


-0.16 


20 


1 


-0.22 


20 


2 


+0.12 


20 


3 


-0.10 


30 


3 


+0.10 


40 


4 


+0.08 


40 


5 


-0.10 



Table 4: Load of power AP^ i (p.u.) for simulation in Scenario 2. +AP^ means a step of required 
power, hence a decrease of the frequency deviation Aw^ and therefore an increase of the power 
reference AP^/^ . 
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1.3 Scenario 3 



We consider the power network described in Scenario 2 and disconnect the area 4, hence obtaining 
the areas connected as in Figure 3. We will simulate Scenario 3 using load steps specified in Table 
5. 
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Figure 3: Power network system of Scenario 3 



Step time 


Area i 




5 


1 


+0.12 


15 


2 


-0.15 


20 


5 


+0.20 


40 


2 


+0.15 


40 


3 


+0.13 


40 


5 


-0.20 



Table 5: Load of power AP^ (p.u.) for simulation in Scenario 3. +APi Ji means a step of required 
power, hence a decrease of the frequency deviation Au^ and therefore an increase of the power 
reference AP re f i . 



2 Design of the AGC layer for a power network using MPC 

The goal of the Benchmark is to design the AGC layer for the scenarios introduced in Section 1. 
Different control schemes will be compared with the centralized MPC scheme described next. For 
a given Scenario, at time t we solve the centralized optimization problem 

F N (x(t)) : (3a) 

t+N-l 

min V (\\x(k)-x \\ Q + \\u(k)-u \\ R ) + \\x(t + N)-x°\\ s ) (3b) 

u(t:t+N-l) 

k—t 

x(k + 1) = Ax(k) + Bu{k) + LAP L (t) ke0:N-l (3c) 

x{k) e X keO: N -I (3d) 

u(k) e U keO: N -I (3c) 

x(N) eX f (3f) 
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and then apply AP re f = u(0). We note that the cost function depend upon x° and u° that are 
defined as xQ = (0, 0, AP^, APj,J and uP, = APj^. The constraints X and U in (3d) and 
(3e) are obtained from constraints listed in Table 2. In the cost function (3b) we set N = 15, 
Q = diag(Qi, . . . , Qm) and R = diag(i?i, . . . , Rm), where 
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and R j = 10. 



Weights Qi and Ri have been chosen in order to penalize the angular displacement A0i and to 
penalize slow reactions to power load steps. Since the power transfer between areas i and j is 
given by 

AP tieiJ (k) = PijiAOiik) - Mj(k)) (4) 

the first requirement also penalizes huge power transfers. 

In order to guarantee the stability of the closed loop system, we design the matrix S and the 
terminal constraint set X/ in three different ways. 

• S is full (M PC full): we compute the symmetric positive-definite matrix S and the static 
state-feedback auxiliary control law K aux x, by maximizing the volume of the ellipsoid de- 
scribed by S inside the state constraints while fulfilling the matrix inequality (A+BK aux )' S(A- 
BK aux ) -S<-Q- K' aux RK aux . 

• S is block diagonal (MPCdiag): we compute the decentralized symmetric positive-definite 
matrix S and the decentralized static state-feedback auxiliary control law K aux x, K aux = 
diag(A'i, . . . , Km) by maximizing the volume of the ellipsoid described by S inside the state 
constraints while fulfilling the matrix inequality (^4 + BK aux )' S(A + BK aux ) — S < —Q — 

K n „,^RI\ aux . 



a ax 



• Zero terminal constraint (M PC zero): we set S = and Xy = 
2.1 Performance criteria 

We propose the following performance criteria for evaluating different control schemes. 

• rj -index 



where T S i m is the time of the simulation. From (5), rj is a weighted average of the error be- 
tween the real state and the equilibrium state and between the real input and the equilibrium 
input. 

• $ -index 

T Bim -i M 

* = ^r- E EE i Ai W*)i T « 

k=0 i=l jeA/i 

where T S i m is the time of the simulation and AP t i eij is the power transfer between areas 
i and j defined in (4). This index gives the average power transferred between areas. In 
particular, if the 77-index is equal for two regulators, the best controller is the one that has 
the lower value of 



G 



3 Control Experiments 



We applied the centralized MPC schemes introduced in the previous section to scenarios 1, 2 and 
3. Furthermore, for each scenario we discrctized the continuos system with both discretization 
schemes D and Dss. At time t we solve the optimization problem (3) and then apply the control 
action to the continuos-time system, keeping the value constant between time t and t + 1. If at 
time t the power load increases or decreases, we assume the controller can use this information 
at time t. This means at time t the controller knows exactly the value of AP^ hence can use 
it. We highlight that violation of this assumption can impact considerably on the index r\. In all 
experiments we use T S i m = 100. In Table 6 and 7 the values of the performance parameters r\ and 
<&, respectively, are reported for each control experiment. 





Scenario 1 


Scenario 2 


Scenario 3 




D 


Dss 


D 


Dss 


D 


Dss 


MPC full 


0.0249 


0.0249 


0.0346 


0.0347 


0.0510 


0.0511 


MPCdiag 


0.0249 


0.0249 


0.0346 


0.0347 


0.0510 


0.0511 


MPCzero 


0.0249 


0.0249 


0.0346 


0.0347 


0.0510 


0.0511 



Tabic 6: Values of the performance parameter rj using different centralized MPC schemes for the 
AGC layer. 





Scenario 1 


Scenario 2 


Scenario 3 




D 


Dss 


D 


Dss 


D 


Dss 


MPCfull 


0.0030 


0.0029 


0.0063 


0.0060 


0.0060 


0.0058 


MPCdiag 


0.0030 


0.0029 


0.0063 


0.0061 


0.0060 


0.0058 


MPCzero 


0.0030 


0.0028 


0.0063 


0.0059 


0.0059 


0.0058 



Table 7: Values of the performance parameter $ using different centralized MPC schemes for the 
AGC layer. 



4 Supporting Matlab files 

We provide the Matlab files for the parameters in Table 2 (parameters. m) and for all control 
experiments. Each file .mat of the control experiments contains 

• the matrices of the continuos linear system {Ac, Be, Cc, Dc, Lc); 

• the matrices of the discrctized linear system [A, B, C, D, L, Ts); 

• parameters of the controller (Q, R, S, N, xO, uO); 

• parameters of the control experiment Tsim and deltaPload, where deltaPload corresponds 
to AP L ; 

• the results of the control experiment x, deltaPref , 77 and $, where deltaPref corresponds 
to AP re f. 
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For each Scenario we included also a Simulink model. In particular, one can load the file .mat 
of a control experiment and simulate the power network system given the power load steps and 
the power reference computed through centralized MPC. 

4.1 Example of simulation 

In the following we illustrate how to use the files .mat and the Simulink models through an example. 
Assume we want to simulate Scenario 2 using the discretization Dss and centralized MPC with 
zero terminal constraint (MPCzero). In the folder of each scenario there are six folders labeled 
as [discretization scheme] _[mpc type]. Hence, we have to use files in folder Dss_MPCzero. In this 
folder we can find the data of the required control experiment as dataSim.mat. The previous 
operations are performed with the Matlab commands: 

cd scenario2 

load Dss_MPCzero /dataSim 

We can simulate different scenarios using the Simulink models present in the folder of each scenario. 
For Scenario 2 we then open the file simulatorPNS-AGC-2.mdl. Start a simulation from Simulink 
will produce the results of the control experiments. These steps are performed with the Matlab 
commands: 

open ( ' simulatorPNS_AGC_2 ' ) 
sim ( ' simulatorPNS_AGC_2 ') 



5 Benchmark exercise 

The aim is to design decentralized/distributed controllers for the scenarios described in Section 1. 

Depending on the control technique adopted either D or Dss discretization schemes can be 
chosen. 

The first goal of a distributed AGC layer is to have performance in terms of ?y similar to 
centralized MPC. Matching also the values of $ can be seen as a secondary objective. 

Alternative control schemes will be also ranked according to the degree of decentralization 
of the design process. Ideally, the controller of each area should be designed independently of 
the others and using information from a limited number of other areas. Decentralized design is 
important in PNS because if an area needs to be isolated or a new area is plugged into the network 
one would like to avoid the redesign the whole AGC layer and rather retune just a limited number 
of local controllers in order to guarantee asymptotic stability and constraints satisfaction for the 
whole network. 
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